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During the last two years there has been much activity and 
discussion at scientific meetings associated with Jit International 
Satellite Cloud Climatology Project (ISCCP). The meteorological 
cormiunity is focusing considerable effort towards the goal of more 
complete understanding of the effects of clouds upon global and 
regional radiation budgets. The First ISCCP Regional Experiment 
(FIRE) workshop defined essential research efforts required to accom- 
plish this goal. An important component of this research plan is the 
modeling and parameterization of stratocumulus cloud fields. 

Consistent with this goal, the present investigation has been 
oriented towards the development of accurate stratocumulus cloud para- 
meterization for use in stratocumulus cloud field studies. One aspect 
of this study has been the expression of the radiative properties of 
broken cloudiness in terms of the better understood horizontally homo- 
geneous plane-parallel values. This type of approach is designed to 
facilitate the utlizatlon of these results by GCM's and climate models. 
In particular, cloud liquid water Is held constant as cloud shape and 
cloud cover is varied. This is important since GCM's predict only 
cloud amount and liquid water content. Only in this way can differ- 
ences in cloud shape be compared and related back to the predicted 
values supplied by GCM's. 

The proximity of neighboring clouds in a cumulus cloud field 
significantly alters cloud field radiance through cloud-cloud photon 
interactions and by shading. Previous studies of cloud field radia- 
tive properties have made one or more of the. following assumptions: 

(1) isotropic radiation from cloud sides, (2) neglect of multiple 
scattering of radiation between clouds, and (3) lack of shading by 
neighboring clouds. In contrast, the present study makes none of 
these approximations, but rather uses the Monte Carlo method to 
simulate cloud field radiative properties exactly. 

The supporting paper, "Stratocunulus Cloud Field Reflected Fluxes: 
The Effect of Cloud Shape" by Welch and Wielicki has been accepted for 
publication in the Journal of the Atmospheric Sciences (November or 
December, 1984). Of particular note is the f act that cloud shape has 
been shown to be a variable as Important as cloud optical density and 
cloud aspect ratio in determining the cloud field radiative properties. 
The size and shapes of gaps between clouds determines such behavior. 
This is because clouds have an anisotropic intensity pattern out the 
cloud sides and because elongated cusp-like regions between clouds 
are efficient at allowing photons to propagate to the ground without 
interaction with neighboring clouds. 

A recent paper by Harshvardhan and Thomas (1984), published in 
the 20 August edition of the Journal of Geophysical Research , presented 
a parameterization scheme for a cl oudTield composed of two-dimensicnal 
bar-ciouds. The enclosed paper by Welch and Wielicki generates the 
three-dimensional version of this technique. 



An analysis of this parameterization scheme shows that It rapidly 
decreases in accuracy for solar zenith angles deviating from e 0 = 60°. 

Thi s method else is limi te d to op ti cel 1y th * v. * , i ^ iu*cua The 
present study has been succesful finding a new parameterization scheme 
that is accurate for solar zenith angles ranging from e 0 = 0° to 82° 
and for all cloud optical depths. The method is the result of a detailed 
examination of the radiation patterns exiting the cloud sides. Final 
calculations are in progress, and a paper describing this parameteriza- 
tion scheme is expected to be submitted for publication by the end of 
the year. 


A Monte Carlo sc.ieme including surface albedo also is close to 
completion. However, this aspect of the study has been downgraded in 
priority in order to complete the above parameterization scheme. Con- 
siderably more effort will be required before the parameterization of 
surface reflection is accomplished. 
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ABSTRACT 


Reflected fluxes are calculated for stratocurr.ulus cloud fields as a 
function of sky cover, cloud aspect ratio, and cloud shape. Cloud 
liquid water volume is held invariant as cloud shape ‘is varied so that 
the results can be utilized more effectively by GCM and climate models. 

Of particular significance is the magnitude of the reflected flux 
differences between broken and plane-parallel cloudiness. On the basis 
of required accuracy in the Earth Radiation Budget Experiment (ERBE) 
program, an order of magnitude value of 10 W m -2 is used to estimate 
"significant differences" between plane-parallel and broken cloudiness. 
This limit is exceeded for cloud covers between 10% and 90%, indicating 
that plane-parallel calculations are not satisfactory at most values of 
cloud cover. The choice of cloud shape also leads to large differences 
in reflected fluxes. These differences may be traced to the anisotropic 
intensity pattern out the cloud sides, to the size and shape of the 
"holes" between clouds, and to variations in cloud area as viewe.! from 
the solar direction. 

An empirical relationship for effective cloud cover is given at 
solar zenith angle of e = 60°. This relationship allows for the rela- 
tively accurate (aF = 10 - 15 W m -2 ) computation of broken cloud field 
reflected fluxes from plane-parallel calculations. Although the present 
parameterization is limited to solar zenith angles near e = 60°, this is 
an indication that further work may lead to reasonably accurate esti- 
mates of broken cloud field radiative properties using modified 
plane-parallel calculations, irrespective of assumed cloud shape. 
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1. Introduction 


S' 


The state-of-the-art of cloud modeling in general circulation models 
(GCMs) is prediction of cloud liquid water content and cloud fraction. 
Parameterizations using relative humidity and/or vertical velocity were 
derived by Kasahara and Washington (1971). More recent improvements are 
typified by Slingo (1980) and Donner et al . (1982) for convective 
clouds, and by Ramanathan et al . (1983) for cirrus. G: /en the 
occurrence of cloud in a GCM grid box and atmospheric layer, the 
radiative properties of clouds are predicted using horizontally homoge- 
neous plane-parallel radiation codes. A number of radiation column 
models of varying sophistication have been developed (Geleyn and 
Hollingsworth, 1979; Fouquart and Bonnel , 1980; Hense et al . , 1982; 
Zdunkowski et a 1 . , 1982; Wiscombe et a 1 . , 1984), some of which attempt 
to include multi-layer broken cloudiness. However, comparisons of the 
radiative output of a large scale atmospheric model with satellite 
measurements (Geleyn ef_ al . , 1982) show a number of weaknesses in 
current radiation schemes, especially for cloud parameterizations. 

Each of the above techniques for modeling cloud radiation is based 
upon plane-parallel assumptions, whereas satellite images show that a 
large proportion, if not the majority, of cloud systems are composed of 
broken cloudiness. Studies by McKee and Cox (1974), Davies (1978), 

Aida (1977), and Harshvardhan (1982) clearly show that broken cloudiness 
is characterized by radiation fields quite different from their plane- 
parallel counterparts. In recognition of these and similar concerns, 
cloud climatology programs such as the International Satellite Cloud 
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Climatology Project (ISCCP) and the First ISCCP Regional Experiment 
(FIRE) are devoted to improving cloud parameter statistics such as cloud 
cover, cloud top height, etc. In addition, so as to make these results 
mere applicable to modeling efforts, these programs also are aimed at 
expressing the radiative properties of broken cloudiness in terms of 
plane-parallel values, adjusted by as yet undefined scale factors 
(Rossow et_ aU_, 1934). 

Early attempts to define the properties of broken cloudiness have 
relied upon the radiative characteristics of individual cloud elements 
(Busygin et al . , 1973; McKee and Cox, 1974; and others). Present efforts 
now have shifted to the modeling and prediction of the radiative proper- 
ties of cloud fields. However, parameteri zations need to be developed 
in a way which facilitates their utilization by GCM's and climate 
models. In the present study, cloud liquid water is held constant as 
cloud shape is varied. Differences in cloud shape then can be compared 
and related back to the predicted values supplied by GCM's. 

Preliminary attempts to model cloud field radiative properties by 
Busygin et al . (1973) and Welch and Zdunkowski (1981a, b) did not include 
cloud-cloud interactions, although Aida (1977) included first order 
effects. Howevei , the preximity of neighboring clouds significantly 
alters cloud field radiance through cloud-cloud interactions and 
shading. Investigations of radiative interactions based upon clouds of 
identical size and shape have been reported by Aida (1976, 1977), Gube 
et al . (1980), Bradley (1981), Harst vardhan (1982), Weinman and 
Harshvardhan (1982), and Claussen (1982). Differences between these 
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studies are attributed primarily to the various approximations used: 

(1) isotropic radiation from cloud sides; (2) neglect of multiple 
scattering of radiation between clcuds; and (3) lack of shading by 
.leighboring clouds. In contrast, the present study eliminates these 
approximations. 

The question is raised as to what degree cloud shape and cloud field 
arrangement are important variables in determining cloud field albedos. 
In particular, for broken cloudiness does albedo tend uniformly towards 
the plane-parallel results as cloud amount is increased? Or, can fields 
of almost overcast clouds have significantly different radiative proper- 
ties from those of the plane-parallel variety? Such questions are 
important both for modeling considerations and to the ISCCP and FIRE 
field programs. The present study seeks to provide a first look at 
these questions by examining the effect of cloud shape upon cloud 
radiative interactions and by determining those values of sky cover for 
which the simpler plane-parallel radiation calculations are adequate and 
those values for which the more complicated broken cloudiness com- 
putations are required. 

2. Method of calculation 

The Monte Carlo method is based upon the Markov chain formulism 
by which the interaction of a photon with cloud droplets is independent 
of all previous interactions. This is a stochastic process for which 
the probabilities associated w’th each scattering event are a function 
only of the present state of a photon packet and not its history 
(Marchuk et al . , 1980). With respect to radiative cloud computations. 
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phot an packets are traced through atmosphere and cloud, accounting for 
both absorption and scattering along the path, until they either reach 
the earth's surface or are returned to space. 
a. Monte Carlo computations 

The Monte Carlo method simulates the photon path; between the (i-l) th 
and the i th scattering event, the photon travels distance Si, as 
determined from the solution: 


Ti - 

(1) 

t i = (*j e j s j)i 

J 

(2) 

Si = (?Sj)i 

J 

(3) 


where Ti is the transmissivity (randomly chosen between 0 and 1), ti is the 
optical path for this event, Pj is the volume extinction coefficient of the 
jth cloud subregion, and Sj is the distance traveled through the j^ 11 sub- 
region. Generally, cloud subregions are partitioned into cubic and rec- 
tangular shapes for ease of computation. Within each subregion or "box," 
cloud attenuation coefficients and phase functions are assumed to be 
uniform. However, these values may vary from box to box. 

After traveling distance Si, scattering initiates a new direction 
of travel. A new azimuthal angle y is randomly assigned with value 
between 0 and 2n, representing rotation of the x-y axes about the axis 
z. A second rotation a of the y-z axes about t;, e axis x' is found 
from the solution to 
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( 4 ) 


p(a') sin a 1 da ' 

-i- = Q(a ) 

f Q P(a') Sin a' da' 

where P(a) is the scattering function; Q(a) Is randomly assigned a 
value between 0 and 1, and then the value of a is determined. 

The photon path is traced until it exits a cloud side. Outside of 
the cloud, the photon is followed until it either reaches the ground, 
enters a neighboring cloud, or is traced to the top of the earth's 
atmosphere. Cloud exit position (x,y,z) and exit angles (9,<t>) for 
each photon are saved for subsequent analysis, 
b. Cloud-cloud interacti ons 

The cloud field is assjmed to be composed of clouds of identical 
diameter (D), shape, and internal optical properties, as discussed in 
Section 3. Spacing between cloud centers (S) is varied to simulate 
cloud cover (N) , where 

N « (0/S) 2 (cubic clouds) (5) 

Table 1 shows values of (S/0) used in this study. 

A photon packet exiting from a cloud side may experience multiple 
scatterings between clouds, entering the sides of neighboring clouds 
before finally reaching the ground or reaching the top of the atmo- 
sphere. In this case, the photon is reentered into the original 
cloud at the equivalent position and angle, as shown by primes in 
Fig. 1 (i.e., periodic boundary conditions). 
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c. Cloud geometry 

Figure 2 shows the cloud geometries assumed in this in/estigation. 

In addition to the cubic shape, the convex geometries of hemispheres, 
cylinders, and capped cylinders (cylindrically shaped bases with 
hemispi.erically shaped tops) are included. 

It should be noted that the cloud convex geometry is defined by the 
appropriate functiors; it is not approximated by smaller rectangular- 
shaped boxes. The intersection of a straight line (i.e., the photon 
path) with the convex function allows the entry and exit points for 
each photon to be precisely located. 

Figure 3 shows an example of a cloud field composed of convex 
shaped clouas of diameter D, and with spacing S between cloud centers. 
Cloud cover N for these clouds is given by 

N = n/4 (D/S) 2 • (6) 

for the array shown in Fig. 3a. Table 1 shows the values of cloud cover 
N as a ^unction of cloud separation (S/D). 

The hexagonal cloud field pattern shown in Fig. 3b occasionally is 
observed in nature. In this case, cloud cover is given by 

N = (tt/1/6) (D/S) 2 (7) 

with values given in Table 2. 
d_. Shadowing 

Mutual shadowing (Fig. 4) occurs for cubic and cyl indrical ly shaped 
clouds for which cloud spacing in the plane of the sun is 
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S/0 < 1 + H/D tan 6 0 


(8) 


The shadow region for the convex shaped clouds Is not as extensive as 
for the cubic cloucs. The point P shown In Fig. 4b, at a distance Z 
from cloud top and at a distance x from cloud center, is shaded by its 
neighbor for 

z tan 9 0 > S - (D 2 - 4x 2 )^ 2 . (9) 

Expressions for cloud shadowing using hemispherical (SPH) and capped 
cylinder (SPHCYL) clouds are more complicated and are not given here. 
Nevertheless, similar conditions apply. 
e. E ffective clou d fraction 

In order to be incorporated into GCM and climate models, cloud 
field radiation budgets need to be expressed in some easily parame- 
terized form. Therefore, the r adiative properties of the cloud field 
are expressed in terms of the better understood plane-parallel values, 
i.e., F C f (N)/F p p(N) , where F c f(N) is cloud field flux at sky cover N 
and Fpp(N) is the place-parallel flux at the same value of cloud cover 

Fpp(N) - N F pp (10G%) . (10) 

The ratio F c f (N)/F pp (N) is equivalent to the ratio N e /N described by 
Sparrow and Cess (1978) and Weinman and Harshvardhan (1982), where N e 
is defined as the equivalent cloud fraction of a plane-pa r al lei cloud 
of the same optical tnickness which produces the same flux as from 
the finite cloud array. Analytical expressions for F c f (N)/F p p(N) 
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appropriate to arrays of cubic and cyl in-irlcally shaped clouds are 
given in Appendix A. 
f_._ Reflected flux 

With respect to using the results (Sec. 3) in GCM and climate 
models, one of the key quantities is the dl f *e r ence in energy reflected 
back to space. The difference between broken cloud and plane-parallel 
reflected fluxes is given by 

aF * F c f (N)-Fpp(N) , (IT) 

expressed in W rr 2 . Using Appendix A, this fiux difference may oe 
expressed as 

AF 3 F c u 0 N [( 1+H/D tan 9o)^c “ Rpp] (12) 

where F 0 is incident solar flux, u 0 ’s the cosine of the solar zenith 
angle 0 o , R pp is planc-pai al lei albedo, and R c is cloud albedo, defined 
as the fraction of photons intercepted by the cloud (top and sides) 
which are reflected back to space. An equivalent expression is 

AF = F 0 u 0 R pp N [Fcf (N)/'F pp (N) - 1] . (13) 

The above results are valid only for non-absorbing regions of the 
solar spectrum. A complete integration across the solar spectrum is 
beyond the scope of the present Investigation. However, a fi^st-order 
correction of Eq. (13), including the effects of absorption and 
Rayleigh scattering, is made below. First, Rayleigh scattering is 
approximated following Lac i s and Hansen (1974). Their parameterization. 
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accurate to 1%, provides estimates of Rayleigh scattering, Rp, of 3.8% 
at solar zenith angle 6 o =0° and 6.6% at 0 o =6O°. 

Second, absorption is estimated from the results of Davis et al . 
(1979), who calculated the distribution of solar absorption in finite 
clouds based upon the partitioning of the solar spectrum into three 
regions. The optical constants, water vapor absorption coefficients and 
energy for each region are given by Welch et aj[ (1976). The tropical 
atmosphere given by McClatchey j?t aj._ (1971) was used, with the top of 
the atmosphere assumed at 15 km. Davis aK found values of absorp- 
tion, A5, above their cloud 630 m thick of 299 W m* 2 at 0 o =O° and 
174 W m -2 at 6 o =60°. These values are nearly independent of cloud 
optical thickness. 

Using the above values, Eq. (13) becomes 

AF = [fo^o(l _ ^R)"^b]RppN[Fcf (N)/Fpp(N) - 13 . (14) 

g. Error analysis 

Using the Gaussian approximation to the binomial distribution of Np 
photons, where P is the probability a photon is reflected back to space 
and (1-f ) is the probability that a photon is transmitted (or 
absorbed), then the stardard deviation for a determination of 
reflectance P is given by 

o p - [Pd-PjNp- 1 ] 172 . (15) 

For the present calculations, cloud field reflectance varies from about 

0.6 to about P = 0.85. For the constant value of N p = 5000 photons 
assumed ir «*ese computations, op ranges from 0.005 to 0.007. 
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The corresponding standard deviation for F C f (N)/F pp (N) is given by 

of = °p R pp A c( 9 o)/ A (16) 

where Ac(0 o )/A is given in Appendix A. At solar zenith angle of 0 o =O°, 
of takes values of 0.006 to 0.009, or less than 1%. 

Similarly, the standard deviation for aF (Sec. 2f) in which clouds 
do not shade their neighbors is 

®AF = (FoUoTN (1+H/D tan 0 O ) o p , (17; 

where (F 0 u 0 )'= F 0 u 0 (l-RR)-Ab is orrected for Rayleigh scattering and 
absorption (Sec. 2f). Note tha„ the error in AF for unshaded clouds is 
oroportional to cloud cover N. For shaded clouds, 

°AF = (^oUo)'N 1 / 2 0p , (18) 

with the error in aF proportional to the square root of cloud cover. 

Now the largest value of o p occurs for P = 0.5, so that 

o p < 0.5 N p -i/2 . (19) 

Using this value, the number of photons necessary to provide a given 
accuracy o A p can be estimated from 

N p = C(F 0 Uo)"°*6 N e /o A p] 2 (20) 

where, for cubes. 
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N(l+n/0 tan 0 O 
Ni / 2 , 


» 


no shading 
shading 


( 21 ) 


Note that Fq. (20) may be expressed as 


N p = (signal/noise) 2 , 


( 22 ) 


where the signal is the plane-parallel reflected flux calculated on 
the basis of cloud fraction N e , and the noise is the uncertainty in aF. 
3. Results 

Stratocumulus cloud fields contain clouds of approximately uniform 
thickness. In the following results, all clouds are assumed to be 1 km 
thick, single scattering albedo is taken as w 0 = 0.999, and the phase 
function is approximated by the Henyey-Greenstein function with asym- 
metry factor g = 0.85. Clouds within the field are assumed to all be of 
equal size and shape and to be arrayed in regular patterns of horizon- 
tally infinite extent. These assumptions are compatible with earlier 
studies by Aida (1977), Claussen (1982), and Weinman and Harshvardhan 
(1982). 

Four different cloud shapes are investigated, as shown in Fig. 2. 

For cubic and cyl indrically shaped cloudr, the volume extinction coef- 
ficient is taken to be 0 e = 49 km -1 . Th’s value is chosen in order to 
compare the present results with previous studies. Cloud diameter is 
varied, keeping cloud thickness fixed, in order to study the effects of 
cloud aspect ratio. 

Plank (1969) found cloud aspect ratio of unity as an average value 
for his aircraft observations of cloud fields, but with H/D ranging 
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between values of about 2 and 1/2. Blackmer and Serebreny (1962) found 
H/D <■ 1, and Bradley (1981) observed H/D ranging between 0.07 and 3.5 
with mean value of 0.73. Finally, Hozumi et al . (1982) found that 
aspect ratio increased with increasing cloud depth, where H < 1 for 
small clouds (D < 1 km), and H/D > 1 for large clouds (D > 3 km). 

In the following results, cloud diameter is varied between D = H/2 and 
D = 2H in order to examine the importance of cloud aspect ratio H/D 
upon cloud field radiative characteristics. 

An important aspect of this investigation is the fact that liquid 
water mass in the cloud field is held constant as cloud snape is varied. 
For example, consider a GCM which predicts liquid water mass m A con- 
densed into a given latitude/longitude grid box and vertical layer. 
Distribution of this liquid water into cloud volume v^d determines the 
average cloud liquid water content w^ = m£/v c ],j. For a given cloud 
fractional coverage in the layer, cloud volume varies with assumed cloud 
shape, causing a variation of liquid water content with cloud shape. 
Assuming a specified cloud droplet distribution, clord extinction is 
linearly proportional to liquid water content. Therefore, comparisons 
of cloud fields of different cloud shapes, normalized to conserve liquid 
water mass in the cloud field, is accomplished by adjusting cloud extinc- 
tion coefficient. These values are increased from <3 e = 49 km -1 to 
B e = 59 km -1 and 0 e = 73 km -1 for capped cylinders and hemispheres, 
respectively. 

a_. Cubic cloud reflected fl ux ratios 

With few exceptions, previous studies of individual cloud and 
cloud field radiative properties have been based upon cubic and 
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rectangular geometries. For the purpose of comparison, this cloud 
geometry is examined first. 

The value of F C f (N)/F pp (N) shown in Fig. 5a for solar angle (0 o =O°) 
varies almost linearly with cloud fraction N. The broken cloud field 
produces cloud radiative flux values which are smaller than the plane- 
parallel values. For N = 1, sky cover is complete and the plane- 
parallel limit is reached. Increasing (decreasing) cloud aspect ratio, 
produced by decreasing (increasing) cloud width, leads to greater 
(smaller) variations from the plane-parallel values. A physical 
explanation of why cloud field fluxes F c f(N) vary with aspect ratio 
for fixed value of sky cover N is deferred to the next section with 
discussion of cyl indrical ly shaped clouds. 

Note that calculations are made at each of the points listed in 
Table 1, as well as ac the shadowing limit. Error bars are small (on 
the order of 1%) and are omitted from the figures, with the curves 
smoothed to within these limits. A complete description of tha error 
analysis is given in Section 2g. 

Figure 6a shows a comparison of the present results with those of 
Aida (1977), Claussen (1982), and Weinman and Harshvardhan (1982) for 
0 O = 0°. The present values and those given by Aida (1977) are nearly 
identical. However, Aida included cloud-cloud interactions only from 
the eight nearest neighbors and allowed only a single photon interaction 
between clouds. In contrast, the present results include all possible 
interactions between clouds in the cloud field. The fact that the 
present values and those given by Aida are nearly identical shows that 
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first order interactions are adequate in simulating cloud-cloud 
interactions, at least at solar zenith angle of 6 0 = 0°. Values of 
F c f (N)/Fpp(N) given by Claussen (1982), based upon an approximation 
which neglected shading, are somewhat larger than the present results. 
Values given by Weinman and Harshvardhan (1982), based upon the 
delta-Eddington approximation and assumed diffuse fluxes from the 
cloud sides are somewhat smaller. 

At solar zenith angle of e 0 = 60°, the stratocumulus field produces 
flux values in excess of those predicted for the plane-parallel case 
(Fig. 5b). The largest variations from plane-parallel occur at the 
shadowing limit. After shadowing occurs, the value of F C f(N) approaches 
the value F pp (N) as sky cover N is increased. F C f(N) > F pp (N) at 
e 0 = 60° since a large portion of incoming radiation is incident upon 
the cloud sides, increasing cloud effective area. Increasing cloud 
aspect ratio produces very large deviations from the plane-parallel 
limit at small values of sky cover. Maximum values of F C f (N)/F pp (W) of 
2.4, 1.7, and 1.4 occur at the shadowing limit for cubes at N = 5%, 13%, 
and 23% for H/D = 2, 1, 1/2, respectively. For sky cover of N > 0.3, 
cloud aspect ratio has diminishing impact upon the value of 
F C f (N)/Fpp(N) ; however, at smaller values of N, aspect ratio is the 
dominant factor determining effective cloud cover. 

No intercomparison of results is shown for 0 O = 60° in Fig. 6a. 
Neither Weinman and Harshvardhan (1982) nor Claussen (1982) provide 
values at this solar zenith angle. Aida (1977) gives a value of 65% for 
his closely packed 9-cloud model. This value is much less than the 85% 
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calculated in the present results using Monte Carlo and various plane- 
parallel methods, indicating that there is an error in his computations. 
Tnerefore, at 0 O = 60°, Aida's results have not been intercompared with 
the present calculations. 
b. Cubic cloud flux di fferences 

For climate models, the key quantity is the difference in radiative 
flux reflected back to space between broken and plane-parallel clouds. 

As discussed in Section 2f, a first order correction for absorption and 
Rayleigh scattering can be made. 

However, while it is beyond the scope of the present investigation 
to directly include absorption for all cloud shapes and aspect ratios, 
a single computation was made for cubic clouds with H/D = 1. This com- 
putation used the same partitioning of the solar spectrum as reported by 
Davis et al . (1979), but with volume extinction coefficient cf B e = 49 
km -1 and cloud thickness of 1 km. Fluxes were normalized to cloud top 
so as to avoid differences due to absorption above the clou<. Resulting 
values of F c f (N)/Fpp(N) are nearly identical to those shown in rig. 5a. 
This is due to the fact that absorption values in plane-parallel and 
finite clouds are nearly identical, in agreement with Davis et al . 
(1979). This result provides further .justification to the procedure 
outlined in Section 2f. 

The differences between broken cloud a n d plane-parallel reflected 
fluxes are computed using Eq. (14) and are given in Fig. 5b. These 
values of AF are calculated using plane-parallel values of cioud reflec- 
tance, Rpp, of 0.79 at 6 0 = 0° and 0.85 at 9 0 = 60°. At e 0 *• 0° and for 
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cloud cover of 5G%, broken cloud fields are darker than their plane- 
parallel counterparts, with differences ranging from 30 W nr 2 for clouds 
with aspect ratio H/D = 1/2 t0 about 60 W m -2 for H/D = 2 (tall clouds). 
In contrast, at e = 60° and for cloud cover of 40 - 50%, broken clouds 
intercept additional energy through their sides and are brighter than 
their plane-parallel counterparts. In this range of cloud cover, flux 
differences, AF, range from about 40 W m -2 for aspect ratio H/D = 1/2 to 
about 60 W nr 2 for H/D = 2. The errors in this analysis are calculated 
from Section 2g using Eqs. (17-18). For N = 0.5, errors are about 
o A p = 5 W nr 2 . The curves in Fig. 5b are smoothed to within these 
limits, and the error bars are omitted. 

The question arises as to what a significant difference is between 
broken and plane-parallel cloudiness. For a region 250 km on a side, 
the Earth Radiation Budget Experiment (ERBE) observational program esti- 
mates a required accuracy of between 2 and 14 W m -2 (Barkstrom and Hall, 
1982) for monthly mean va 1 ;S of reflected solar fluxes. Using 10 W m -2 
as an order of magnitude for significant differences in reflected 
fluxes. Fig. 5b suggests that finite cloud and plane-parallel results 
are significantly different for cloud covers ranging between N = 0.1 and 
0.9. While the current study examines these differences for two solar 
zenith angles, a more complete treatment will require averages over sun 
angles throughout the day. The present study points out the large 
potential for such differences. 

Figure 6b gives the comparison of aF with other studies at e 0 = 0°. 
Calculated from the results of Weinman and Harshvardhan (1982), the 
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maximum difference in fluxes, AF, is 30 W m -2 , compared to 40 W nr 2 
calculated from Aida (1977) and the present study; Claussen's (1982) 
calculations result in slightly larger flux differences, 

£. Cylindrically s haped clouds 

For deep, growing clouds which constituted more than half of his 
sample, Bradley (1981) found cloud shape to be nearly cylindrical. 
Figure 7 shows a LANDSAT image of a stratocumulus field with high cloud 
cover and with close cloud spacings. In many cases shown here, the 
clouds may be approximated by circular cross-sectional areas. For the 
sake of intercomparison with cubic clouds, the cloud distribution shown 
in Fig. 3 is used for most of the results presented in this section. 
Nevertheless, Fig. 7 suggests that the hexagonal cloud array is more 
representative for real cloud fields. Results for both the linear 
and hexagonal arrays are given below. 

Figure 8a shows values of F C f (N)/Fpp(N) as a function of cloud 
cover. There are many similarities with the results for cubes shown in 
Fig. 5, but also notible differences. In particular, for the linear 
a r ray at closest cloud spacing (N = 0.79), the curves do not tend 
towards the plane-parallel limit F c f (1)/F pp (l) = 1. In order to 
reach the plane-parallel limt as sky cover is increased, the cylinders 
must be progressively deformed to fill the holes between the clouds. 

In order to gain insight concerning the shape of the curves shown 
in Figs. 5a and 8a, refer to cloud area ratio Aq(60°)/Ac( 0°) as plotted 
in Fig. A. 1 . Note that the shading limit for cylinders occurs for 
N = 0.24, 0.11, and 0.06 for cloud aspect ratios of H/D = 1/2, 1, and 
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2, respectively. At small values of sky cover, before shading occurs, 
the area ratio has tne constant value (1+H/D tan 0 O ) for cubes, but 
the larger values of (1+4 * H/D tan 0 O ) for cylinders. Examination of 
Figs. 1-3, along with simple geometrical analysis, reveals that the 
curved sides of cylinders produce this larger area ratio. After the 
shadowing limit is exceeded, the value of area ratio Aq(60°)/Aq( 0°) 
decreases rapidly with increasing value of sky cover. For large values 
of N, the area ratio varies as N _1 . 

The increase In values of F C f (N)/Fpp(N) at small values of N, up 
to the shadowing limit, is due to cloud-cloud interactions. However, 
the decrease in values after the shadowing limit is reached demonstrates 
the important result that geometrical shadowing effects are dominant 
over increases due to cloud interactions. Weinman and Harshvardhan 
(1982) reached a similar conclusion for studies of two-d'.mensional bar 
clouds. 

Figure 8b shows the corresponding values of aF. At 0 O = 60°, 

F C f(N) exceeds the plane-parallel value Fp p (N) by maximum values of 
75 W nr 2 at sky covers of N = 0.4 - 0.5. These values are about 
15 - 30 W m -2 larger than corresponding values found for cubes 
(Fig. 5b). Choice of aspect ratio and cloud array patterns have a rela- 
tively small (lu-15 W m -2 ) impact upon the values of AF at 0 O = 60°. 
Cloud fields composed of cylinders a r e brighter than those composed of 
cubes because the cylinders have larger effective areas, as shown in 
Fig. .M. 

The situation is diametrical’y opposite at 6 0 = 0°, for which case 
aspect ratio is an extremely important variablf. Maximum value of AF 
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Increases from 45 W nr ; ' for H/D = 1/2 to about 90 W m -2 for H/D » 1 to 
about 130 W m -2 for H/D = 2. These differences are much greater than 
those predicted for cloud fields composed of cubic geometries. Therefore, 
for 8 0 = 0° and for constant value of cloud cover, a cloud field com- 
posed of cylinders is considerably darker than one composed of cubes. 
Differences between cylinders and cubes are smallest for flat clouds 
(H/D = 1/2) and greatest for tall clouds (H/D = 2). 
d. "Holes" between clouds 

A discussion of the mechanisms producing the large differences 
between cubes and cylinders follows. Consider the case of N = 0.7 for 
which these differences are large. In this instance, cubes have separa- 
tion distance between cloud centers of S/D = 1.2; that is, for cloua 
diameter of 1 km, there are gaps between the clouds of 0.2 km (Fig. 9). 

For the same cloud cover and cloud diameter, cylinders have a separation 
distance of S/D = 1.05; that is, the minimum gap between clouds is 0.05 
km. Now consider the diagonal distance between the clouds shown in 
Fig. 9. For cubes, this distance is 0.28 km and for cylinders, 

0.48 km, or about twice as far. 

First, for aspect ratio H/D = 1, zenith angle e 0 = 0° and cloud 
cover N = 0.7, cubes and cylinders produce nearly identical zenith angle 
scattering patterns for photons initially exiting the cloud (Fig. 10A, 
First Exit). These are the scattering patterns associated with non- 
interacting clouds. About 33%-34% of incident energy is scattered out 
the cloud tops and about 2% is transmitted out the cloud bases. Of the 
65% of incident energy scattered out the cloud sides, about half exits 
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In the upper quarter of the cloud; only about 2% exits from the lower 
quarter of the clouc 1 sides. Of the total energy out the cloud sides, 
about 25% is scattered into the upper hemisphere, and about 39%-40% 
scattered into the lower hemisphere. At zenith angle of e 0 = 60° 

(Fig. 10B, First Exit), cylinders scatter a smaller proportion of Inci- 
dent energy out the clojd tops than do cubes; however, cylinders scatter 
more energy than cubes out the cloud sides. This is due to the fact 
that cyMnders have a smaller average path length through the cloud at 
e 0 > 0° than do cubes. 

For 0 O = 0°, Fig. iCA (Final Exit) shows the photon scattering 
pattern on final exit from the cloud. This is the energy distribution 
after all cloud-cloud interactions have taken place. Note that a large 
redistribution of energy occurs through the cloud-cloud interactions. 

The percentage of incident energy scattered out the cloud sides drops 
from about 65% for photons first exiting the cloud to about 26% for 
protons on final exit from cubic clouds. Cylinders scatter a much 
larger percentage, about 35%-36%, out their sides. The result is that 
cubes scatter more energy out their tops and bases than do cylinders. 

A similar condition applies at zenith angles of e 0 = 60° (Fig. 10B, 

Final Exit). 

Consideration of these differences between cubes and cylinders 
provides an explanation of how the shape of holes between clouds ys 
an important role in cloud reflected fluxes. First, the scattering pat- 
tern on “First Exit" shows that photons travel more frequently in the 
downward than upward direction. Second, the majority of photons exiting 
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from the cloud sides leave from the upper quarter. Finally, the angular 
distribution pattern for photons out the cloud sides Is anisotropic, 
with few photons exiting at zenith angles of C°-30° and 150°-180°. 

Perhaps most striking is the small number of photons on final exit 
from the cubes in the angular bin 90°-120°, and the large number of pho- 
tons exiting in the angular bin 30°-60°. This is due to the fact that 
the majority of photons exit from the upper quarter of tne cloud side. 
Those photons traveling in the 30°-60° angular bin more easily avoid 
further cloud-cloud interactions, whereas photons in the 120°-150° angu- 
lar bin almost never avoid entering a neighboring cloud. In contrast, 
the relatively large gaps between clouds for the cylinders lead both to 
greater scattering from the clouds' sides and to a mure symmetrical 
distribution of scattered photons. However, in all cases the scattering 
pattern remains highly anisotropic, with few photons in the 0°-30° and 
150°-180° angular bins. 

These angular patterns demonstrate that photons exiting cloud 
sides require considerable distance between neighbors in order to escape 
interaction with them. Since the greatest distance between clouds at 
N * 0.7 is about 0.3 km for cubes, and about 0.5 km for cylinders, 
many more photons are able to penetrate the cloud gaps without inter- 
action from cylinders than from cubes. The result Is that cylinders 
reflect less total incident energy for equal values of sky cover than 
do cubes. 

As aspect ratio increases, cloud width decreases for clouds of 
equal vertical extent; cloud separation distance then decreases for 
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constant values of sky cover N. If the angular cxltance pattern were 
Invariant with aspect ratio, then Increased Interactions would make 
clouds of larger aspect ratio brighter. However, this Is not the case. 
Increased numbers of photons exit the cloud sides as aspect ratio 
Increases, and exit angles become more strongly peaked the downward 
direction. These considerations more than compensate for the decrees^ 
in cloud separation, thereby reducing cloud reflectance for octn cubes 
and cylinders as aspect ratio Increases. 

This analysis shows that the size of gaps between clouds strongly 
Influences the total energy reflected by a cloud field at constant sky 
cover N. Choica of the hexagonal cloud pattern, rather than the 
linear array used above, dec; eases the average gap d. stance between 
cylindrically shaped clouds. The result is a brighter cloud field 
for the same cloud cover. 

Close examination of Fig. 7 shows a number of examples of the 
hexagonal cloud array. Gaps between clouds also tend to be cusp-shaped, 
indicative of convex shaped clouas; linear patterns would indicate 
clouds with straight sides. As cloud cover N is increased towards 
N * 1 , the clouds must become progressively deformed in order for 
F c f(l) to approach the plane-parallel limit Fpp(l). Examination of the 
LANDSAT image In Fig. 7 indicates that clouds seldom are exactly cir- 
cular in shape and almost never are cubic. Therefore, results between 
those for cubes and cylinders may be more appropriate. 
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e. Clouds with hemispherical tops 

While most clouds are neither cylindrical nor cubic in shape, 
they also seldom have the flat cloud tops assumed in the previous 
results. The purpose of this section is to examine the effect of 
clouds having tops with the hemispherical shape. Aspect ratio is held 
constant with a value of H/D = 1 for these cases. 

At 9 o = 0°, the value of F c f (N)/Fpp(N) for hemispheres (Fig. 11a) 
is similar to that for cylinders with aspect ratio H/D * 1/2. For 
capped cylinders (H/D = 1), values of F c f (N)/Fp p (N) increase more 
rapidly than do cylinders (H/D = 1). However, due to the much smaller 
values of effective area (Fig. Al) at 0 O = 60°, capped cylinders and 
hemispheres are significantly darker than their cylindrical or cubic 
counterparts at the same value of cloud cover N. Differences are 
largest at the small values of cloud cover. Shading limits occur at 
N = 0.14 for capped v./i'inders and N = 0.35 for hemispheres 

At 9 0 = 0°, capped cylinders have values of AF (Fig. lib) which are 
smaller than for cylinders, but larger than for cubes. These values of 
aF are nearly identical to those of cylinders with the small aspect 
ratio H/D = 1/2. 

At 0 O = 60°, the presence of curved cloud tops tends to minimize 
differences between plane-parallel and broken cloudiness. 
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4. Discussion 

Enhanced values of reflected flux at solar zenith angle of 60° 
are largely due to geometrical effects. As shown by the ratio 
N e /N = F C f (N)/Fpp(N) , effective cloud cover N e becomes much larger than 
actual cloud cover N. An attempt to correct for these geometrical 
effects can be made by redefining AF as 

AF' = F cf (N) - F p p(N') , (23) 

where N' is an adjusted cloud cover. Taking N' = Ac(6)/A as defined in 
Appendix A, AF' becomes 

AF* = [F C f(N)/F pp (N) - N'/N] NR pp [F 0 u 0 (l-R R )-A b ] (24) 

Figure 12 shows the values of AF' for solar zenith angle of 60°. 
Comparison of Figs. 5b and 8b with Fig. 12 shows that the values of AF' 
are negative and as large in magnitude as the values of aF. This 
suggests that an average between aF and aF' may be more accurate. 
Therefore, the empirical relationship 

AF" = [F cf (N)/F pp (N) - (l+N7N)/2] NR pp [F 0 u 0 (l-R R )-A b ] (25) 

is applied and is shown in Fig. 11. This kind of empirical fit also 
has been suggested by Harshvardhan and Thomas (1984). For cubes and 
cylinders, the value of AF" is negative at small values of cloud cover 
and positive at larger values of N. However, the magnitude of aF" is 
less than 15 W nr 2 for cylinders and less than 10 W m -2 for cubes. For 
capped cylinders and hemispheres, the value of AF remains negative for 
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all values of cloud cover. The magnitude of AF is less than 10 W m -2 
for capped cylinders and less than about 17 W m -2 for hemispheres. 
Therefore, the empirical fitting procedure allows for a relatively 
accurate method of obtaining cloud field reflected fluxes from plane- 
parallel calculations at solar zenith angle of 60°. 

A caution is in order, lest the reader believe that a reliable 
method has been developed to calculate cloud field reflected fluxes. 

The accuracy of the AF" fit is limited to solar zenith angles near 60° 
and has not been established at cloud optical depths other than t = 49. 
Note that at solar zenith angle of 0", N' = N, so that no improvement in 
the results is possible using the empirical method. As solar zenith 
angle is decreased from 60°, the error in aF" using the empirical fit 
must steadily increase (not shown in the figures). 

5. Conclusions 

The present study has been careful to preserve cloud liquid water 
volume as cloud shape was varied. It is assumed that cloud optical 
volume scales linearly with liquid water volume. In this way, the 
results give the differences in solar reflected fluxes which would be 
predicted for a cloud field with specified cloud fraction and cloud 
liquid water volume but varying cloud geometry. In particular, cloud 
fields with plane-parallel, cubic, cylindrical, hemispheric, and 
hemisphericiil-capped-cylinder geometries are compared for a range of 
cloud cover. 

Values of F C f (N)/F p p(N) shown in the figures give the scale factors 
which relate the albedos of broken cloudiness to the more familiar 
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plane-parallel values. At 0 O = 0°, the albedos of broken cloud fields 
are distinctly smaller than corresponding plane-parallel values; at 
e 0 - 60°, the broken cloud field albedos are larger, witn maximum 
differences occurring at the shadowing lim^.. 

However, perhaps of greater significance is the magnitude of the 
reflected flux differences between broken and plane-parallel cloudiness. 
For fixed solar zenith angle, the value of 10 W nr 2 is used here as an 
estimate of "significant differences," chosen on the basis of detect- 
ability by broad-band flux measurements. This limit of aF < 10 W m~ 2 
occurs for N < 0.1 and for N > 0.9, indicating that plane-parallel 
calculations a^e not satisfactory at most values of sky cover. 

For the larger values of sky cover, differences in AF can become 
very large for some cloud shapes. This is due to the fact that the 
scale factor F C f (N)/F pp (N) is relatively "flat" as a function of N 
at e 0 = o° for convex-shaped clouds and does not approach the planer- 
parallel limit as N * 1. Deformation from this circular cross section 
occurs in real clouds (Fig. 7), indicating that cloud albedos in 
between those predicted for cubes and cylinders may be more appro- 
priate. The presence of convex cloud tops also serves to decrease the 
difference between plane-parallel and broken clouds and between cubes 
and cylinders. The net result is that cloud computations based upon 
cubic and rectangular geometries may provide reasonable estimations of 
cloud field albedo, at least in selected cases, even though real clouds 
do not have these geometries. Such agreement also is suggested by the 
observations of McKee et al . (1983). 


In any case, the present results show that cloud shape and the 
size of the gap*, between clouds are important variables. Clouds have an 
anisotropic intensity pattern out their sides and elongated cusp-like 
regions between clouds are efficient at allowing photons to propagate to 
the ground without interaction with neighboring clouds. Therefore, the 
presence of small holes, as opposed to long thin lines, between clouds 
may have a large impact upon cloud albedo. At large values of sky 
cover, these gaps may tend to act as "light pipes" in increasing cloud 
transmissivity and decreasing cloud albedo, producing results signifi- 
cantly different from those predicted by plane-parallel models. 

At solar zenith angle of 60°, geometrical effects dominate, since 
effective cloud cover becomes much larger than actual cloud cover. An 
empirical relationship for effective cloud fraction was found to produce 
errors in aF less than 10 W m -2 for cubes and capped cylinders and less 
than about 15 W nr 2 for cylinders and hemispheres. This technique, 
based upon a similar procedure by Harshvardhan and Thomas (1984), allows 
for the relatively accurate computation of broken cloud field reflected 
fluxes from plane-parallel calculations. However, the use of this 
empirical relationship is limited, since as solar zenith angle decreases 
from 60°, the error in aF using this relationship increases. At solar 
zenith angle of 0°, effective loud fraction becomes cloud cover N. 
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APPENDIX A 


The value N e is defined as the eaui valent cloud fraction of a 
plane-parallel cloud of the same vertical optical thickness (horizontally 
averaged over the cloud field) which produces the same flux as from the 
finite cloud array, 

F C f « N e F pp (Al) 

where the plane-parallel flux is expressed as 

Fpp = Rpp Fo • (A2) 

Plane-parallel albedo is R p p„ F 0 is incident solar flux, and u 0 is the 
cosine of the solar zenith angle e 0 . Flux from the cloud field array is 
given by 

F c f(N) « R c f Uo F 0 , (A3) 

where cloud field albedo R c f is 

Rcf = «c ^C(N,9 0 )/A , (A4) 

taking ground albedo as zero. Cloud albedo R c is the fraction of photons 
intercepted by the cloud (top and sides) that are reflected back to space, 
A c( N » 9 o) the projection of cloud area onto a horizontal surface as 

viewed from the solar direction, and A~S 2 is the total area of a unit cell 
surrounding the cloud. 
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Cubic clouds. For cubic clouds of thickness H and width D, cloud 


intercepted c ea is 


Ac(N,0 o ) 


D 2 + DH tan e 0 , [S/D > 1 + H/D tan e 0 ] 


DS, 


[S/D < 1 + H/D tan e 0 , shading]. 


(A5) 


Combining Eqs. (A4) and (A5), cloud field albedo becomes 


Ref « 


R c N (1 + H/D tan e 0 ), [S/D > 1 + H/D tan e 0 ] 

R c N l/2 , [S/D < 1 + H/D tan e 0 , shading]. 


(A6) 


Scaling the plane-parallel flux (Eq. A2) by cloud cover N produces 


Fpp(N) = N Fp p . (A7) 

Finally, combining Eqs. (A1-A7), we obtain the albedo enhancement factor 
(or effective cloud fraction ratio): 

1(1 + H/D tan 9 0 ) , [no shading] 

Fcf(N)/F pp (N) = N e /N = R c Rpn { ( A8 ) 

L N- 1 / 2 , [shading]. 

This is the amount by which cloud field reflectance is increased above 
or decreased below the equivalent planp-paral lei value, Fpp(N). 

Cylindrical clouds . For cyl indrical ly shaped clouds, unshaded area is 
given by 

Ac(N,e 0 ) ■ "D 2 /4 + DH tan e 0 , [S/D > 1 + H/D tan e 0 ] (A9) 

Partially shaded area is given by 
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Ac(Mo) 


OS, for [S/0 < H/D tan 0 O ] 

OH tan e 0 + 0 2 /0 sin** ( S / 0 . h/D tan 0 O ) 

+ (S - H tan 0 o )/2[O 2 - (S - H tan 0 o ) 2 ] 1/2 , 
[H/D tan 0 O < S/0 < 1 + H/D tan 0 O ] 


(A10) 


Combining the above expressions leads to the albedo enhancement factor for 
cylindrically shaped clouds: 


F c f(N)/F pp (N) = RcRjJi 


(1 + 4 /it H/D tan e 0 ), /or [S/D > 1 + H/D tan 0 o ] 
4/tt[H/D tan 0 O + 1/2 sin" 1 (S/D - H/D tan 0 O ) 

(S/D - H/D tan 0 O )/2(1-(S/D - H/D tan 0 O ) 2 ) 1/2 ], 


| for [H/D tan 0 O < S/D < i + H/D tan 0 O ] 
^ 4/tt S/D for [S/D < H/D tan e 0 ] 


(All) 


Capped cylinders and hemispheres . For capped cylinders (SPHCYL) and 
for hemispherical 1y shaped clouds (SPH), the albedo enhancement factor is 
more mathematically complicated. Th:refore, it is given here in the form 

F C f (N)Fpp(N) * RcR-^A C (N,0o)[Ac(N,O o )] (A12) 

where area enhancement factors, Ac(N,60°)/Aq(N, 0°) , are shown in Fig. A1 as 
a function of N for each cloud shape. 
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APPENDIX B 


The results presented previously express cloud field radiative 
fluxes in terms of plane parallel values. These computations are pro- 
portioned to cloud field area enhancement values (e.g.. Fig. Al) and to 
the ratio of interacting to non-interacting cloud field reflectance 
values. Let values R and Rj represent cloud reflectance values for 
interacting and non-interacting clouds, respectively; then ratio R/Rj 
is the reflectance enhancement factor representing the effect of pho- 
tons, exiting from one cloud, which enter a neighboring cloud. 
Reflectance enhancement ratios are given in Fig. fil for cloud fields 
composed of cubes, cylinders, and hemi spherical ly topped clouds, 
respectively. 

Reflectance ratios for cubic clouds at solar zenith angles of 
e 0 = 0° and 60° are shown in Fig. Bl(a,b). These values increase with 
increasing value of cloud cover. For complete cloud cover (N = 1), 
cloud interactions increase cloud albedo 27%, ab^j that which would be 
found for completely r jn-interacting clouds. These values are nearly 
identical to values reported by Aida (1977) and Weinman and Harshvardhan 
(1982). For aspect ratios of H/D = 2 and 1/2, maximum increase in 
cloud albedo is 57% and 12%, respectively, over that produced by non- 
interacting cloud fields. Clearly, cloud aspect ratio Is a variable 
of great importance to cloud field albedo. 

At e 0 ■ ^°, cloud interactions enhance cloud albedo by 44%, 67%, 
and 89% for !I/D * 1/2, 1, and 2, respectively. The discontinuity in 
values shown in Fig. Bl(b) occu r s at the cloud mutual shading limit. 
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Figure Bl(c,d) shows corresponding values for cylinders. Albedo 
enhancement ratio at 0 O = 0° is considerably smaller for all aspect 
ratios than corresponding values of cubes. Differences are due to the 
fact that photons existing from cloud sides have less chance of inter- 
acting with cylindrically shaped clouds than with cubic shaped clouds. 
Maximum sky cover is N = 0.79 for the checkerboard array and N = 0.91 
for the hexagonal array. Note that only the single case of aspect 
ratio unity is shown for the hexagonal array (HEX). At 9 0 = 60° > 
values of R/R L for cylinders are similar to those values for cubes, 
indicating again tnat shadowing is a dominant influence on values of 
cloud field albedo. 

Figure Bl(e,f) shows similar values of R/P.j for hemispheres and 
capped cylinders. Cloud extinction coefficients are adjusted from 
0 e = 49 km -1 to values of B e = 73 km -1 and 59 km -1 for hemispheres 
and capped cylinders, respectively, so as to preserve cloud optical 
volume. Capped cylinders have larger values of R/R t at 0 O = 0° than 
do the flat-topped cylinders. In contrast, hemispheres produce values 
of R/Rj similar to those for cylirderr with aspect ratio H/P = 1/2. 

At 0 O = 60°, both capped cylinders and hemispheres have values of 
R/Rj significantly smaller than those produced by cylinders. These 
differences are associated with shading effects produced by the curved 
upper surfaces. 
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FIGURE CAPTIONS 


Fig. 1: Cloud field array composed of cubic shaped clouds of diameter 

(width) D with cloud spacing S. Paths of three photons are 
shown which interact with neighboring clouds; primes indicate 
the reinsertion points for these photons with the original 
cloud. 

Fig. 2: The four cloud geometries used in this investigation are 

shown here. 

Fig. 3; Cloud field arrays for convex shaped clouds of diameter D 
with cloud spacing S, arrayed in A) linear and B) hexagonal 
patterns. 

Fig. 4: Shading by neighbors for A) cubic shaped clouds; and 

B) cyl indrical ly shaped clouds. 

Fig. 5: A) Ratio of cubic cloud field radiative fluxes to plane- 

parallel fluxes as a function of cloud cover; B) The dif- 
ference between cubic cloud field and plane-parallel 
radiative fluxes in W m -2 as a function of cloud cover. 

Fig. 6: Comparison of present and previous results at 9 0 = 0° for 

cloud field radiative fluxes. A) Ratio of cubic cloud field 
radiative fluxes to plane-parallel fluxes as a function of 
cloud cover; B) The difference between cubic cloud field and 
plane-parallel radiative fluxes as a function of cloud cover. 

Fig. 7: LANOSAT Multispectral Scanner (MSS) image with 80 m 

resolution taken off the coast of southern California on 
23 September 1979. Horizontal scale is 185 km and vertical 
scale is 170 km. 
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Fig. 8: A) Ratio of cylindrical cloud field radiative fluxes scaled to 

plane-parallel fluxes as a function of cloud cover; B) The 
difference between cylindrical cloud field and plane-parallel 
radiative fluxes i.i W m -2 as a function of cloud cover. 

"HEX" refers to cylinders (H/D = 1) arrayed in the hexagonal 
pattern (Fig. 3b). 

Fig. 9: Sizes (in km) of the gaps between clouds (at N = 0.7) for 

cloud fields composed of cubes and cylinders. 

Fig. 10: The angular distribution of radiation out the cloud sides 
at N = 0.7 for cubes and cylinders at solar zenith angles 
of A) 9 0 = 0° and B) 60°. The angular pattern is divided 
into segments each 30° in zenith angle width. The concentric 
circles represent exiting energy of 5%, 10%, 15%, 20%, and 
25% of total incident energy for photons first exiting the 
cloud, and 2% 4%, 6%, and 8% of total incident energy for 
photons finally exiting the cloud (including cloud-cloud 
interactions). The boxes to the left of each half-circle 
show the percentage of total incident energy which exits 
each quarter of the cloud sides, from top to bottom. 

Fig. 11: A) Ratio of cloud field radiative fluxes scaled to plane- 
parallel fluxes as a function of cloud cover for capped 
cylinders and hemispheres; B) The difference between cloud 
field and plane-parallel radiative fluxes as a function of 
cloud cover. 
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The difference between cloud field and plane-parallel 
radiative fluxes as a function of cloud cover calculated 
on the basis of the A) AF'; and B) AF" approximations. 

Area enhancement ratios as a function of cloud cover N for 
the four cloud shapes shown in Fig. 2. Cloud aspect ratio is 
H/D = 1. "HEX" tefe.s to cylindrically shaped clouds arrayed 
in the hexagonal pattern (Fig. 3b). 

Reflectance enhanceicnt factors as a function of sky cover 
at solar zenith angles e 0 = 0° and 60° for (a,b) cubes; 

(c,d) cylifK- s; (e,f) capped cylinders and hemispheres. 
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S/0 

CUBIC CL0U0S 

CONVEX-SHAPED CLOUDS 

CONVEX-SHAPED CLOUDS 

1.00 

100.0 

Square Array 
78.5 
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90.7 

1.05 

91.0 
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82.2 
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